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Recent experiments have demonstrated that it is possible to alter the dispersion of a medium 
without significantly altering its absorption or refractive index and that this may be done while a 
wave propagates through the medium. This possibility opens up a new set of potential experiments 
to the field of nonstationary optics. We consider the basic kinetics of waves propagating through 
a medium whose group and phase velocities are a function of position and time. We compare the 
dynamics of waves propagating through two homogeneous media, one a nondispersive medium with 
a time dependent phase velocity, and one a dispersive medium with a time dependent group velocity 
and show that new dynamic effects accompany new kinetic ones. 
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I. INTRODUCTION 

The terms nonstationary, inhomogeneous, and disper- 
sive describe media whose properties vary with time, po- 
sition, and frequency, respectively. In this paper, we con- 
sider transformations wrought upon waves as they prop- 
agate through nonstationary, inhomogeneous media with 
variable dispersion. Our interest in this topic is moti- 
vated by experimental advances in the manipulation of 
dispersion through nonlinear optics. Recent experiments 
have demonstrated that it is possible to achieve either 
strong normal dispersion (in the case "slow-light media 
[HIS]" ) or strong anomalous dispersion (in the case of "su- 
per luminal media [3]") in nearly transparent media us- 
ing nonlinear optical effects. In either case, the strength 
of the dispersion, as experienced by a probe beam, de- 
pends on the intensity of one or more controllable aux- 
iliary beams. If auxiliary intensities are changed while 
the probe beam is in transit, then the dispersion, as per- 
ceived by the probe beam, is time dependent. These 
media, which we will refer to collectively as controllably 
dispersive media, thus form a newly accessible class of 
nonstationary media. Although the nonstationarity of a 
slow light medium has been used experimentally to adi- 
abatically transform a "slow light" pulse to a "stopped" 
one [2JI3], an explicit connection between this application 
and the field of nonstationary electromagnetics has not, 
to our knowledge, been made. 

The field of nonstationary electromagnetics, although 
much less developed than its inhomogeneous counterpart, 
is now a venerable one. In 1958, Morgenthaler intro- 
duced propagation equations for electromagnetic waves 
in an isotropic nondispersive medium whose permittivity 
and permeability were allowed to vary with time (but not 
with space) |S] . Although the propagation equations are 
not generally analytically soluble, he found useful solu- 
tions assuming, on the one hand, a step functional time 
dependence and, on the other hand, an adiabatic time 
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dependence. He found, among other things, that the fre- 
quency of a wave would vary with the permittivity and 
permeability in such a way that wavelength would be 
preserved. 

Over the next two decades, much of the progress of the 
field was in the Soviet Union, with a particular focus on 
the behavior of light in plasmas. We refer the reader to 
Stepanov and colleagues (SI [7] for an introduction to this 
literature and highlight here two aspects of it. First, the 
two simplified cases (step function and adiabatic changes) 
used by Morgenthaler continued to be used extensively, 
and were generalized from a purely nonstationary pic- 
ture to one that allowed inhomogeneities to propagate 
at a fixed velocity according to the traveling wave law 
(/ = f{x~vt)) OIH]. This picture allows a single frame- 
work to unite the purely nonstationary case {v oo) 
P7] , a 1-dimensional purely inhomogeneous case [v — Q), 
and the case of an inhomogeneity (for example, an ion- 
ization front) moving at any velocity in between these 
two extremes. Second, a distinction between "kinetic" 
and "dynamic" phenomena was found to be useful in dif- 
ferentiating between those results of the theory that are 
general to all linear wave phenomena and those that are 
specific to particular media O |H1 12] • We employ both of 
these concepts in this paper. 

The advent of the laser, and particularly the pulsed 
laser, introduced a new way to dynamically alter the 
characteristics of a medium. Although they were appar- 
ently unaware of previous work in nonstationary electro- 
magnetics, a few authors noted the possibility of using 
light to alter the phase velocity of propagating electro- 
magnetic waves. In 1977, Lampe, Ott, and Walker noted 
that an ionizing laser might be swept across a gas to cre- 
ate a superluminal ionization front that could interact 
with a microwave pulse in a nonstationary fashion |10) . 
In 1988, Wilks, Dawson, and Mori examined the prob- 
lem of a wave propagating through a medium which is 
then quickly ionized by a high-intensity ultrashort pulse 
[llj . Large frequency shifts were soon realized for waves 
reflecting off of a relativistically propagating ionization 
front [12]. These works founded a sub-literature sur- 
rounding nonstationary effects in plasmas [TMIS] . We 
point the reader to Shvartsburg, who reviewed some of 
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the recent work on nonstationary effects in plasmas [T7] . 
We also highlight the work of Avitzour and Shvets, whose 
2008 paper relies on controlled dispersion in nonstation- 
ary media. They proposed a method for altering the 
spectral width of a pulse without changing its central 
frequency using controlled dispersion in a nonstationary 
magnetized plasma [T5]. We will look at this method 
from a kinetic perspective in the next section. 

In parallel with the more recent work on nonstationary 
plasmas, several authors have explored electromagnetic 
propagation through nonstationary media abstractly, in 
the tradition of Morgenthaler [THHH] . Dodonov, Klimov, 
and Nikonov studied the quantization of a linear, non- 
stationary medium and were able to quantitatively re- 
late temporal changes in dielectric permittivity to pho- 
ton generation . Biancalana, Amann, Uskov, and Or- 
eilly treated a 1-dimensional nondispersive nonstationary 
medium by introducing a transmission matrix for mov- 
ing interfaces. Using this matrix they generalized Bragg 
reflection to propagating interfaces and showed that tem- 
poral periodicity leads to k-vector bandgaps just as spa- 
tial periodicity leads to frequency bandgaps [20j. Budko 
has found an analogy between a nonstationary nondis- 
persive medium and an expanding Universe [21] . 

One surprising fact about recent works, both those fo- 
cusing on nonstationary plasmas piMTB] and those deal- 
ing with abstract, nondispersive, linear media [T^HS], is 
that they make no reference to the substantial Soviet lit- 
erature from the 1960s and 1970s. One consequence of 
the disjointedness of which this fact is symptomatic is 
that basic principles have had to be rediscovered, mul- 
tiple times, in incidental and typically incomplete ways. 
For example, a basic principle of wave propagation in 
nonstationary linear media is the importance and gener- 
ality of kinetic effects. Simple spatial and temporal sym- 
metries play an important role in determining the effects 
of nonstationary propagating waves. Aspects of these 
symmetries have been rediscovered and used in many 
contexts HI HH dHl UHl 121] ■ 

We hope that this paper, in addition to providing a 
simple theoretical basis for treating a new class of non- 
stationary dispersive media, will demonstrate the utility 
of explicitly recognizing kinetic constraints and of main- 
taining a clear distinction between kinetic and dynamic 
aspects of wave behavior. Section|Tl]is dedicated to an ex- 
ploration of the interplay between dispersion, nonstation- 
arity, and inhomogeneity from a strictly kinetic perspec- 
tive. We introduce three basic kinetic relationships rep- 
resenting three basic preserved symmetries. Along with 
appropriate dispersion relationships, they allow for the 
derivation of all the subsequent kinetic results. Although 
these relationships have not, to our knowledge, been ex- 
plicitly presented together before, they have been implicit 
in many previously derived results. We show that they 
lead to Snell's Law for a motionless interface and to ap- 
propriate Doppler shifts for a moving one, that they lead 
to the preservation of wavelength for homogeneous non- 
stationary media and to the preservation of frequency in 




FIG. 1: An incident plane wave with k-vector ki interacts 
with a planar interface. This results in a refracted wave with 
k-vector kr and a transmitted wave with k-vector kt- For a 
non-dispersive and isotropic medium, the refractive index for 
the reflected wave (rir) is the same as that for the incident 
wave (rii). 



inhomogeneous stationary media, how they give rise to 
Biancalana et aZ.'s "generalized frequency" and give in- 
sight to Avitzour and Shvets' proposal for compressing 
pulse spectrum without altering carrier frequency. We 
also find some effects that are described here for what 
may be the first time. For example, we find that disper- 
sion modulates the frequency response of a wave to tem- 
poral changes in the refractive index and that Doppler 
reflections may lead to large changes in pulse bandwidth 
when the moving interface approaches the group veloc- 
ity. We emphasize that although we have derived these 
results with controllably dispersive nonlinear optical me- 
dia in mind, the results apply to any type of wave prop- 
agating through any linear medium. 

Unlike kinetic effects, the dynamic aspects of wave be- 
havior depend on the specific microscopic details of in- 
teractions between nonstationary media and fields. In 
Section [mi we develop boundary conditions appropriate 
to an idealized version of the nonlinear optical media cur- 
rently used to achieve exotic dispersion. Unsurprisingly, 
we find that these boundary conditions are fundamen- 
tally different from those posited for nondispersive media. 
We compare pulse transformations wrought by a chang- 
ing phase velocity in a specific, homogeneous, nondisper- 
sive medium (originally described by Morgenthaler ^), 
with those wrought by a change in group velocity in our 
idealized dispersive medium. We note many interesting 
similarities and differences in their effects on twenty key 
quantities such as energy, momentum, and photon den- 
sity. 
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II. KINETICS OF DISPERSIVE, 
INHOMOGENEOUS, NONSTATIONARY MEDIA 

Consider a plane wave incident on an infinite planar 
interface, as drawn in Fig. [l] Its wave vector and fre- 
quency describe periodic translational symmetries in 3 
dimensions of space and 1 of time. A planar interface 
breaks symmetry in only one dimension, i.e., the one par- 
allel to its normal. If an interface is motionless, then its 
normal is purely spatial; in the geometry of Fig. [l] we 
could write n — x. The normal of a moving interface in 
3-1-1 dimensions is partially rotated into time; if we intro- 
duce the temporal unit vector ct and a temporal rotation 
angle 0, then we may write h = cos(f>x + sin (pet. The 
angle (j) is related to the velocity v and the normalized 
velocity f3 by 



c 



tan ( 



(1) 



All the other translational symmetries of the incident 
plane wave are preserved under reflection and transmis- 
sion. Thus, 



ky2 ^yr ^yt-i 



kzi ^zr ^zt7 



and 



c c 



(2) 
(3) 



(4) 



where the geometry is again defined by Fig. [T] These 
three equations, combined with appropriate dispersion 
relationships, provide a basis for the kinetic aspects of 
wave behavior. 

Equations (2]), ([3|, and ^ represent quantities that 
are unchanged by a planar interface. Extending this con- 
cept, we see that if these quantities are not altered by a 
single interface, they cannot be altered by multiple sim- 
ilar interfaces ( "similar" means here that they share the 
same velocity and the same normal) . By progressively in- 
creasing the number and variety of similar interfaces, we 
can similarly see that the same quantities are preserved 
under interaction with any pattern of medium parame- 
ters that obeys the traveling wave law. 



(5) 



When the planar disturbance described by the func- 
tion / in Eq. ([5| is more complicated than a single inter- 
face, the description of the wave response in terms of a 
single incident, reflected, and transmitted wave becomes 
inadequate. In order to express conserved quantities ap- 
propriately for this more general case, we now introduce 
the convention that a line appearing over a given quan- 
tity means that the quantity is invariant under planar 
changes in the material parameters. Using this conven- 
tion, we extend Equations (|2]), ([3]), and Q to 

ky, (6) 



and 



cos(0) — kx sin(0). 



(7) 



(8) 



Despite their lack of an explicit equality, we refer these 
statements equations in acknowledgment of the fact that 
they are, in fact, compressed equalities. 



A. Snell's Law and Doppler shifts 

For a motionless interface, /3 = and Eq. Q reduces 
to the statement that frequency is preserved: 



(9) 



Combining this with Eq. ^ and the expansion kz 
nujsm9/c yields 



smf 



sm f 



sm Ut, 



(10) 



which is Snell's Law (note that Snell's Law applies to the 
reflected wave as well as the transmitted wave). 

If we allow the interface to move at a velocity v, fre- 
quency is no longer conserved, and Snell's Law will no 
longer hold [5S] . If both sides of the interface are disper- 
sionless, then a simple extension of Snell's Law obtained 
by combining Eq. ^ and Eq. ^ may be useful: 



rit sm Of 



n, sin 9i 



1 — rit cos 9t[3 1— Ui cos 9if3 



(11) 



One interesting fact about this formula is that it shows 
that there are some cases where there two sets of angles 
and frequencies for the transmitted wave [29] . When dis- 
persion is important, Eq. (11) loses its utility because 9t 



and Tit both become frequency dependent. We may still 
find angle, refractive index, and frequency using both re- 
lationships 



nt{uJt)i^t sin9t ^ ni{uji)uji sin( 



(12) 



and 



ujt{l - nt{ujt) cos 9tl3) ^ uJi{l - ni{LL>i) cos 9^/3), (13) 

combined with an explicit function for nt{uit) to calculate 
the properties of the transmitted wave. We will explore 
the effects of dispersion in a simplified case, that of the 
transmitted wave at normal incidence, in the next sub- 
section. 

The relationships for the reflected wave. 



nr{uJt)ujr sin9r — ni{uji)uji sin^^ 



(14) 



and 



Wr(l — nr{LOr) COS 0^/3) = Wi(l — ni{uji) COS 9ij3), (15) 
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combined with an explicitly function for nr{cOr) appear 
identical to those for the transmitted wave. However, the 
behavior is different in this case because the medium that 
defines is the same as the medium that defines rii and 
because the sign of cos 9 is different in the refiective case. 
Numerical investigation suggests that the Doppler effect 



we obtain 



grows large when 



V for an approaching wall. In 



Subsection |II E| we will show analytically that this is the 
case for a wave that is normal to the interface. For now, 
we note that for normal incidence on a moving interface 
in a vacuum, cosfl^ = — 1, cos0, — 1, and n„ — rij — 1. 



Then Eq. ( 15 ) gives the standard reflective Doppler shift. 



1-13 
1 + /3' 



(16) 



Note that in the geometry of Fig. [T] a positive value for 
/3 means that the interface is receding. 



B. Transmission at normal incidence 

In controllably dispersive media, large changes in 
group velocity may be accompanied by minute changes 
in phase velocity. In addition, the phase velocity changes 
are gradual in the sense that they occur over several car- 
rier frequency cycles. Under these conditions, reflection 
is negligible. At normal incidence, cos^i = cos 6*4 = 1. 
Then Eq. ([s]) simplifles to 



— cos (j) — ksinc 
c 



(17) 



where we have removed the subscript x from in ac- 
knowledgment of the fact that at normal incidence k — 
kx- This mixture of wave vector and frequency is the con- 
served quantity that Biancalana et al. found and referred 
to as the generalized frequency [201 . When i/i = 0, the 
medium is stationary and frequency is preserved. When 
<j) = ±7r/2, the medium is homogeneous and wave vec- 
tor is preserved. In between these extremes, a mixture 
between wave vector and frequency is preserved. 

The preservation of the generalized frequency is ac- 
companied by a second constraint: 



k 



(18) 



where Vp is the phase velocity. When the phase velocity 
changes, the ratio Lo/k must change as well. Together, 
Eqs. (17) and (18) determine the frequency and wave 



vector responses. We now look at the impact of disper- 
sion upon these responses. Consider a medium for which 
the spectral dependence and the time dependence of the 
refractive index are independent so that we may apply 
n = rioix, t) + {5n/ {Suj)){uj — uJo) for frequencies near Wq. 
From 



d 
drir, 



COS( 



fcsin 



in0^ 



0, 



(19) 



du 

dUn 



uj sm<p 
cos d — Jin sin < 



(20) 



where Ug = n + ujSn/Suj. Since tan0 = v/c, we can also 
write 



du) 

dun 



ujv/c 

1 - V/Vg ' 



(21) 



which emphasizes the large frequency shifts expected 
when the interface velocity v approaches the group ve- 
locity Vg. 



C. Preservation of a generalized bandwidth in 1 + 1 
dimensions 

The group velocity is defined using Vg = duj/dk. When 
group velocity dispersion may be neglected, either be- 
cause it is small, or because we are limiting ourselves to 
a narrow spectrum, the approximation 



(22) 



becomes useful. As formalized in Eq. (22), when a nar 



row band pulse experiences a change in group velocity, 
either its spectral bandwidth Aw, or its spatial band- 
width Ak, or both, must change. Extending Eq. (17 1, 
we find a generalized bandwidth. 



Auj 



Ak sin ( 



(23) 



that is preserved when a pulse interacts with an inter- 
face moving at the velocity defined by </> and Eq ([I]). 
Eq. ( 23 ) shows that the division of change between spec- 



tral and spatial bandwidth depends on their comparative 
magnitudes and on the interface velocity. 

The preserved generalized bandwidth provides us with 
a simple tool that we can use to analyze one particular 
aspect of the stopped light experiment of Liu, Dutton, 
Behroozi, and Hau [2]. In the Liu experiment, a coupling 
beam was used to control the group index perceived by 
a probe beam in real time. However, this coupling beam 
had a finite speed and was copropagating with the probe 
beam. Because of this, changes to the group index per- 
ceived by the probe beam were not instantaneous, but 
moved at nearly the speed of light. In their report [2], 
and in earlier theoretical work [22^ , the speed of the cou- 
pling beam was treated as effectively infinite, because it 
was so much faster than the slowed advancement of the 
probe pulse. Using the preserved generalized bandwidth, 
we now show that treating the velocity of the coupling 
beam as essentially infinite was justified. 

When the pulse described by Liu et al. entered the slow 
light medium, it crossed a spatial interface, and (j) was 
0. The preserved quantity associated with this case was 
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FIG. 2: Field energy density verses spatial extent of a pulse. 
S: when a pulse undergoes a spatial change in group index 
(ug) it scales longitudinally but not temporally. T: when a 
pulse undergoes a temporal change in (rig), it changes tem- 
porally but not longitudinally. By mixing spatial and tempo- 
ral interfaces, it is possible to scale the pulse bandwidth and 
longitudinal extent without otherwise changing the envelope 
shape. 



Auj/c, meaning that spectral bandwidth was preserved. 
The spatial bandwidth (Afc) and the group index (rv) 
were increased by 7 orders of magnitude, and Eq. ( 22 ) 
remained satisfied. 



While the new foreshortened pulse was propagating 
through the slow-light medium, the coupling beam was 
then turned off. The speed of the coupling beam was ap- 
proximately c, so that (f) was approximately 7r/4. Thus, 
the preserved generalized frequency was approximately 

{Auj/c — Ak)/V'2., which appears to distribute changes 
evenly between spectral and spatial bandwidths. How- 
ever, the slow speed of the probe envelope was reflected 
by the fact that Ak was larger than Aujjc by a factor 
of 10''. According to the preservation of the general- 
ized bandwidth, Aw/c would have gone to zero while Afc 
changed by 1 part in 10^. This is the size of the error 
introduced by assuming that the coupling speed was ef- 
fectively infinite. 

The end result of the Liu experiment was a reconsti- 
tution of the original pulse. In the next section we use 
the concept of preserved generalized bandwidth to show 
how the pulse may be spectrally compressed or expanded 
using transmission in controllably dispersive media. 



D. Pulse compression/decompression using 
controlled dispersion 

The generalized bandwidth is conserved across any 
number of interfaces so long as they have the same veloc- 
ity. By mixing interfaces of different velocities, we can 
change that bandwidth. In Fig. [2j a pulse begins in a 



medium with a group index rig of 1, transitions to one 
with rig = 2, and then back to = 1. When both tran- 
sitions are spatial (S corresponds to a motionless planar 
interface) or when both are temporal (T corresponds to 
time-dependent but spatially independent change to the 
medium-the infinite velocity case) the pulse ends up in its 
original form. By mixing spatial and temporal transfor- 
mations, the pulse can either be compressed or extended. 
In this way the pulse bandwidth may be compressed or 
decompressed while the pulse envelope is scaled longitu- 
dinally but otherwise unchanged. 

This mechanism for pulse compression and decompres- 
sion was recently proposed and demonstrated in simula- 
tion for a plasma geometries, where the group velocity 
may be controlled dynamically via an applied magnetic 
field [TBI . We note here that the Liu experiment demon- 
strates all the capabilities necessary to perform pulse 
compression or expansion using nonlinear optics [5]. 



E. Dispersive modulation of the 1-D Doppler shift 

In a dispersive medium, Doppler shifts are difficult to 
predict intuitively because of the interdependence be- 
tween refractive index and frequency. Here we find an 
intuitive formula for the change in frequency in the case 
where the group velocity differs from the phase veloc- 
ity, but where the range of frequencies concerning us is 
sufficiently narrow that we may neglect group velocity 
dispersion. Expanding Eq. Q for the case of reflection 
at normal incidence, we get 



(24) 



To derive a Doppler shift, we must introduce an explicit 
dispersion relationship n — n{uj). In general, the result- 
ing equation is analytically intractable though a numer- 
ical solution is readily found. An interesting argument 
for the physical relevance of the concept of group veloc- 
ity to Doppler shifts in dispersive media is that if we re- 
quire that the group velocity be constant with frequency 
a simple analytical solution results. Requiring that n^, 
the group index, remain constant with frequency gives a 



refractive index of the form n 



-b/uj^ where h can be 



any real constant. Applying this dispersion relationship 



to Eq. (24) yields the exact expression. 



A 



-2n- 



(25) 



where /3 = u/c, A = uj^ — uji, ujr + i^i ~ 2ti;, and n = 
n{uj). Eq. (25) is attractive because of its symmetry but 



incomplete because it depends on the quantity w, which 
it does not resolve. Rewriting Eq. (25) in terms oiuji and 
A 



A —2n(LiJi)v I c 



(26) 
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which clearly reveals the role of dispersion in determining 
the expected Doppler shift. Eq. (26 1 has the interesting 



interpretation that the Doppler shift will be large when 
the interface is approaching (if Vg > 0) or receding (if 



< 0) with a speed 



This result is exact when 



the group velocity is constant with frequency. 



Dispersion and the 1-D Doppler effect on pulse 
bandwidth 



Intuitively, we expect that a pulse reflecting from an in- 
terface whose velocity is close to that of the pulse will be 
temporally stretched if the interface is receding and com- 
pressed if the interface is approaching. We expect this 
temporal stretching/compressing should have a corollary 
decrease or increase of the spectral bandwidth, Aw. We 
now show that this is the case. 



Adapting expression (231, we find 



AkjV — Ao 



AkrV, 



(27) 



where Aw^, Afc^, and Awr, Afc^ are the spectral and spa- 
tial bandwidths of the incident and reflected pulses, and 
V is the interface velocity. Assuming that the initial and 
final pulses are sufficiently narrow such that each may be 
associated with a specific group velocity, we may combine 



CVo 



Au/Ak with Eq. (271 to obtain 



AoJr 

Auji 



1 - 



1 



(28) 



where Vgi is the group velocity of the incident pulse and 
Vgr is the group velocity of the receding pulse. As always, 
a positive sign for v means a receding interface and a 
positive sign for Vg means a positive group velocity. Thus, 
as expected, when a pulse reflects from an interface that 
is receding at roughly the pulse group velocity of the 
approaching pulse, the reflected pulse has a narrowed 
bandwidth. When a pulse reflects from a boundary that 
approaches at roughly the group velocity of the reflected 
pulse, the reflected pulse is temporally compressed and 
has a broad bandwidth. 



G. How group velocity modulates the frequency 
response of a pulse to a time-dependent refractive 
index 

When interface velocity, v, is zero, (j) is also zero and lj 
is conserved. When velocity is infinite, </) = ±7r/2. The 
refractive index is a function of time only and k is con- 
served. In this way, the interface model is easily extended 
to treat the case of a time-dependent (but spatially homo- 
geneous) medium. Since k is conserved, frequency must 
compensate for a change in refractive index according to 



(29) 



where we take the medium to have an initial refractive 
index Uo and the pulse to have an initial central frequency 
LJo- If n were driven to 0, uj would be up-converted by 
many orders of magnitude. The changing refractive index 
associated with plasma generation has been proposed for 
use in frequency up-conversion in several works |11H131 
[T51 [TCI [231 [23]. On a more pedestrian level, frequency 
response to a time-dependent refractive index is also at 
the heart of the function of acousto-optic and electro- 
optic modulators. 

Interestingly, the frequency response to a changing re- 
fractive index can be modulated by the group velocity 
of a pulse in a medium. To see this, we begin with a 
dispersive medium with a nondispersive time-dependence 
whose refractive index takes the form n = ni){uj)+t5n/ 5t. 
(This might be taken to model, for example, a gas that is 
a mixture of two atoms, one with resonance that is near 
Lo, the time-dependent frequency of the wave traveling 
through the medium, and another whose resonances are 
all spectrally distant but whose concentration changes 
with time.) Then, from Eq. (29) we find 



did 

~dt 



Lo Sn 

Tin St ' 



(30) 



where n„ is the group refractive index, defined by n„ = 



n + uj6n/Suj. One interpretation of Eq. (30) is that the 
frequency response of a pulse to a temporal change in 
refractive index is proportional to the distance the pulse 
travels while the change occurs. 



III. A COMPARISON OF TWO 
NONSTATIONARY HOMOGENEOUS MEDIA: 
KINETICS AND DYNAMICS 

In this section we seek to augment our understand- 
ing of the effects of a time-dependent group velocity in 
a dispersive medium by comparing them with the effects 
of a time-dependent phase velocity in a no77.dispersive 
medium. To do so, we now define two scenarios, which 
we call Scenarios A and B. In both scenarios, we imag- 
ine a spectrally narrow pulse traveling through a time- 
dependent medium and ask how that pulse is transformed 
as the medium changes with time. 

In Scenario A, we take the medium to be nondispersive, 
spatially homogeneous, and isotropic. We assume that 
changes in the medium are adiabatic, spatially homoge- 
neous, and isotropic. We allow e and fi to vary gradually 
and independently over time. We assume that loss at all 
times is negligible, meaning that e and fi are taken to be 
real for the narrow spectral range of the pulse. Because 
the medium is nondispersive, the group velocity is equal 
to the phase velocity. 

In Scenario B, we take the medium to be dispersive, 
spatially homogeneous, and isotropic. We assume that 
changes in the medium are adiabatic, spatially homo- 
geneous, and isotropic. We fix e and /i for the central 
frequency, but allow their slopes to vary gradually and 
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Number Quantity 


Symbol 


Morgenthaler 


Liu 


1 


Phase velocity 


Vp 


Vpf/Vpi 


1 


2 


Group velocity 




Vpf/Vpi 


Vgf/Vg, 


3 


Permittivity 


t 




1 


4 


Permeability 






1 


5 


Wave number 


k 


1 


1 


6 


Spatial bandwidth 


Ak 


1 


1 


7 


Central frequency 




Vpf/Vpi 


1 


8 


Spectral bandwidth 


Auj 


Vpf/Vpi 




9 


Total energy density 


Ut 


Vpf/Vpi 


1 


10 


Electric displacement 


D 


1 


-J^gshgi 


11 


Magnetic induction 


B 


1 




12 


Electric field 


E 




^/Vgf/Vgi 


13 


Magnetizing field 


H 




\/'Vgfhgi 


14 


Photon energy 


hu 


Vpf/Vpi 


1 


15 


Photon number density 


N 


1 


1 


16 


Poynting vector 


S 


{Vpf/Vpif 


'"g^hg^ 


17 


Abraham momentum density 


E X H/c^ 


{Vpf/Vpif 


'"gfhg^ 


18 


Minkowski momentum density D x B 


1 


'"gShg^ 


19 


Canonical momentum density 


N{nfiw/c) 


1 


1 


20 


Angular momentum density 


N{h) 


1 


1 



TABLE I: A comparison of narrow-pulse transformations by two dynamic media. 



independently over time. We take loss at all times to be 
negligible, meaning that e and /i are taken to be real for 
the narrow spectral range of the pulse. 

A. Kinetics: basic symmetries 

The spatial homogeneity of Scenarios A and B implies 
that 

kf = ±fc„ (31) 

and the adiabaticity of the two scenarios eliminates re- 
flection, requiring that we use only the positive sign. Spa- 
tial homogeneity also implies that 

Afc/ = Afc, (32) 

for the pulses of the two scenarios. Since |fc| = nuj/c and 
Afc « nglS.u)/c, changes in the phase and group velocity 
imply proportional changes in the central frequency and 
pulse bandwidth respectively. Since Ug = n + LoSn/Sio, 
changes in n imply changes in Ug for a nondispersive 
medium and therefore changes in n lead to changes in 
bandwidth as well. 

To continue our comparison, we now specify particu- 
lar boundary conditions for the two changing media. In 
doing so, we lose the generality of Section [llj but gain 
the ability to examine the effects of changes to specific, 
time-dependent media on basic quantities like energy and 
momentum densities and reflection and refraction coeffi- 
cients. 



We take our two media from previous works. The non- 
dispersive medium is taken from Morgenthaler [S]. The 
dispersive medium is an idealization of the dynamic EIT 
medium used in the Liu experiment [2]- 



B. Dynamics: boundary conditions and reflection 

1. Morgenthaler' s nondispersive medium 

Morgenthaler allowed for the electric permittivity (e) 
and magnetic permeability (/^) of the medium to be func- 
tions of time, but not of space or direction or frequency. 
That is, e(r, k, a;,t) = e(i), and ^(r, k, a;,i) = ^{t). He 
assumed that charge and flux, and therefore the elec- 
tric displacement (D) and the magnetic induction (B), 
would be conserved at an interface boundary. Under this 
assumption, the reflection and transmission coefficients 
for the electric field are given by [S] 

and 

|E,| 2 {er ^ y firer J ' ^ ' 

In any Morgenthaler medium, the magnitude of the mag- 
netic field is given by |H| = |E|/ry, where 7] = -\//i/e is 
the impedance. In general, there is a reflected wave. 
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FIG. 3: EIT energy diagram. The coupling beam and the 
medium may be taken together to form an efTective medium 
as felt by the probe beam. The properties of the effective 
medium depend on the state of the coupling beam. 

We consider the case where changes to the medium are 
adiabatic and homogeneous. In this case the reflected 
power vanishes. Using the fact that what is preserved 
across one interface must be conserved across many par- 
allel interfaces (see the beginning of Section |ll]) , we find 
that we may take the boundary conditions across an adi- 
abatic change to be 

D/=D, (35) 

and 

B/ = B,. (36) 

The values of E and H at any time can then be found 
via e and /i. 

2. Liu's dispersive medium 

The Liu experiment involves a pulse interacting with 
an EIT medium across a narrow spectral range around 
the EIT resonance. The system as a whole comprises two 
beams of light interacting with a cloud of cold atoms. 
The two beams are tuned to two connected atomic tran- 
sitions which share an excited state. The stronger beam, 
referred to as the coupling beam in Figure [3| is consid- 
ered in combination with the cloud of atoms to form the 
EIT medium, as seen by the weaker beam, which we call 
the probe beam. When EIT is established, the atoms are 
in a coherent dark state, so-called because it does not 
couple to the probe beam. The dark state may be rep- 
resented as a combination of electronic levels |1) and |2) 
as 

_ l^c|l) ~ ^p|2) exp[i{kp - kc) -f- i{ujp - ujc)t] 
^^l + 

(37) 

If the probe beam is sufficiently weak such that the 
atomic number density remains significantly larger than 
the probe photon number density, then the EIT medium 
will be approximately linear in probe intensity, so that 
the only sizable nonlinear effects of the EIT medium are 



the effects of the coupling beam on the probe beam. Sur- 
prisingly, this linear regime may be maintained even as 
the coupling beam is completely attenuated \F]\ . 

For our nondispersive medium, we idealize the Liu ex- 
periment [2j by neglecting decoherence and other losses 
and require that total pulse energy be conserved under 
changes to the group velocity. Thus, 

utf=uti, (38) 

where Uf represents the total energy density associated 
with the wave including not only that portion that is 
in field form but also that portion that is coherently 
stored in the medium and is a cycle-averaged quantity. 
For a narrow-band pulse propagating through a lossless 
but dispersive medium, we may define a non-dispersive 
energy density (u„) which is related to the dispersive en- 
ergy (ut) by 

n 

Un = Ut — . (39) 

Ug 

In terms of the macroscopic fields, u„ is given, for the 
isotropic medium we are considering, by 

Un^l {eE^ + k^H^) , (40) 

where E'^ and are time averaged quantities. (Note 
that Un matches the standard form for the electromag- 
netic energy density in a lossless, dispersionless, macro- 
scopic medium |25] . When dispersion becomes negligible, 
Ug ss n, and Ut reduces to this form.) 

If, as we have assumed for our comparison, changes 
in the group velocity occur around a constant refractive 
index, then we may take the quantity to be fixed. 

We will also assume that the impedance at line center is 
unchanged, so that the ratio between H and E remains 
fixed. Because the refractive index at line center will also 
remain fixed, all macroscopic electromagnetic fields scale 
together. Therefore, all terms quadratic in any pair of 
fields scale with the group velocity. The direction of the 
fields is unchanged by the homogeneous change in the 
isotropic medium. We summarize this as 

FfV^ = FV^, (41) 

where F could be any of the four macroscopic electro- 
magnetic fields. 



C. Transformations wrought under changing phase 
and group velocities in the 'Morgenthaler' and 'Liu' 
media 

We are now able to perform a more complete compar- 
ison between scenarios A and B. In Table |l) we list 20 
different quantities and show how they are scaled as the 
media of Scenarios A (Morgenthaler) and B (Liu) change. 
Quantities 1-4 summarize the difference between the two 
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scenarios, emphasizing that the phase velocity and group 
velocity change in Scenario A while only the group veloc- 
ity changes in Scenario B. The scaling of quantities 5-8 
depends only on symmetry (see Section [ll]) and therefore 
are general consequences of the scaling of quantities 1-4. 
The scaling of quantities 9-13 may be taken as expres- 
sions of the boundary conditions for the two media. 

The remaining quantities (14-20) can be found in sim- 
ple ways from a knowledge of the scaling of the first 13 
quantities. The per photon energy (14) scales with fre- 
quency, which scales as in line 7. That the photon num- 
ber density (15) is conserved is related to our assumption 
of no loss. The Poynting vector (quantity 16) scales as 
the square of the phase velocity for the non-dispersive 
case. This quadratic dependence combines the linear de- 
pendence of the total energy density on the phase ve- 
locity with the fact that the group velocity is equal to 
the phase velocity for a nondispersive medium. Because 
the dispersive case assumes that the phase velocity re- 
mains unchanged, the Poynting vector of the Liu medium 
feels only the effect of the change in group velocity. The 
Minkowski momentum (18) is conserved in the nondis- 
persive case, but is proportional to the group velocity for 
the nondispersive case. This second dependence under- 
scores the fact that the Minkowski momentum, as defined 
by D X B, does not give the total momentum when dis- 
persion is taken into account. The canonical momentum 
(19), which docs give the total momentum, is conserved 
for both media. Because the total photon number is con- 
served for both media, so is the total angular momentum 
density associated with the pulse for both media. We 
note that although the total angular momentum does not 
change, its distribution between field and medium does. 

In concluding this section we emphasize three points. 
First, to understand the effects of a particular medium 
on a pulse, it is necessary to specify appropriate bound- 
ary conditions that depend on the details of the way in 
which the changing medium interacts with the field. Sec- 
ond, boundary conditions that are reasonable for a non- 
dispersive medium (Scenario A) may fail for a simple 
dispersive medium (Scenario B) for the particular reason 
that changes in dispersion imply a shifting of energy from 
field to medium. It then becomes necessary to model dif- 
ferent changes on a case-by-case basis. Finally, a general 
comparison of the scaling of different quantities for Sce- 
nario A and Scenario B shows that the changes wrought 
in a controUably dispersive nonstationary medium are 
qualitatively different from those wrought in a nonsta- 
tionary nondispersive medium. In this sense, controllably 
dispersive media open up a new and potentially fruitful 
niche in nonstationary electromagnetics. 

IV. CONCLUSION 

Controllably dispersive media are an experimental re- 
ality |2}i4J . All the requirements necessary to perform in- 



teresting nonstationary experiments in controllably dis- 
persive media have been demonstrated for one experi- 
mental system [21 S], are near at hand for another [3], 
and appear to be feasible for others [T51 [55] . In this pa- 
per we have outlined a basic theory for the treatment 
of controllably dispersive media from the perspective of 
nonstationary electromagnetics. 

Although allowing the properties of a medium to de- 
pend not only on position but also on both time and fre- 
quency opens up a large parameter space, the behavior of 
waves in this space is subject to basic kinetic constraints 
imposed by symmetry. We have explicitly defined the re- 
quirements of symmetry for the case of a plane wave in- 
teracting with a moving planar interface. We have justi- 
fied these requirements by rederiving well-established re- 
sults, such as Snell's Law and the free space Doppler shift, 
and used them to simply derive analytical descriptions 
of more esoteric phenomena, such as the effect of disper- 
sion on the reflective Doppler shift. Using the outlined 
kinetic theory, we have provided an analysis that cor- 
roborates a recent proposal [T8| for the compression and 
decompression of pulse envelopes in magnetized plasmas 
using strictly nonstationary effects and have shown that 
the proposed mechanism may be applied in any control- 
lably dispersive medium. We have introduced a simple 
expression that shows how a pulse may be significantly 
compressed or decompressed under Doppler reflection- 
even if the Doppler shifts themselves are small. Finally, 
we have shown that dispersion modulates the frequency 
response of a wave to a temporal change in the refractive 
index-the response is damped for slow light media, and 
amplified for fast light media. 

Unlike kinetic effects, which rely only on simple sym- 
metries, dynamic effects are determined by boundary 
conditions which depend upon particular details of the 
medium at hand. We have compared dynamic effects 
in two simple, homogeneous, adiabatically nonstationary 
media-one which is nonstationary in phase velocity and 
a second which is nonstationary in group velocity. We 
compared the effective boundary conditions associated 
with these media and the changes wrought in spectrally 
narrow pulses propagating through them. We have found 
substantial differences in both the boundary conditions 
and the effect of the nonstationary medium on funda- 
mental quantities such as energy, and electromagnetic 
momentum. These differences suggest that it may be 
useful to consider controllably dispersive media as a new 
and potentially fruitful category of nonstationary elec- 
tromagnetic media. 
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